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PREFACE 


The  work  described  in  this  report  was  authorized  under  the  In-House  Laboratory 
Independent  Research  grant  from  the  U.S.  Army  Edgewood  Chemical  Biological  Center.  This 
work  was  started  in  January  2008  and  completed  in  March  2011. 

The  use  of  either  trade  or  manufacturers’  names  in  this  report  does  not  constitute 
an  official  endorsement  of  any  commercial  products.  This  report  may  not  be  cited  for  purposes 
of  advertisement. 

This  report  has  been  approved  for  public  release. 
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REGRESSION  MODEL  FOR  MODTRAN  WITH  APPLICATIONS 
TO  INACTIVATION  OF  MICROBES  SUSPENDED  IN  THE  ATMOSPHERE 
BY  SOLAR  ULTRAVIOLET  RADIATION 


1.  INTRODUCTION 

The  inactivation  of  airborne  viral  pathogens  is  relevant  to  public  health  and 
biodefense.  In  a  recent  paper  (/),  we  described  MODerate  resolution  atmospheric 
TRANsmission  (MODTRAN),  a  software  program  developed  by  Spectral  Sciences,  Inc. 
(Burlington,  MA)  and  the  U.S.  Air  Force  Research  Laboratory  (Wright-Patterson  Air  Force 
Base,  OH)  (2).  Using  MODTRAN,  we  explored  the  effects  of  altitude,  sea-level  meteorological 
visibility  distance,  and  surface  (ground)  albedo  on  the  inactivation  of  viruses  due  to  direct  and 
diffuse  solar  UV  in  the  wavelength  range  280-320  nm  (known  as  UVB)  in  a  clcar-sky 
atmosphere  (2). 

Virucidal  solar  radiation  consists  of  a  direct  solar  UV  and  a  diffuse  UV 
component,  which  is  produced  by  the  multiple  scattering  of  light  by  ambient  atmospheric 
aerosols  and  gases  and  the  multiple  reflections  of  UV  by  the  underlying  surface  or  ground.  In 
this  study,  wc  extend  the  regression  model  results,  which  were  limited  to  diffuse  and  total 
irradiances  in  our  previous  work  (/),  and  we  outline  the  general  procedure  to  compute  the 
regression  coefficients.  The  procedure  may  be  useful  for  setting  other  MODTRAN  irradiance 
computations  in  the  framework  of  a  simple-to-use  regression  model. 

In  our  previous  paper  (/),  we  discussed  the  difference  between  plane  (horizontal) 
irradiances  and  spherical  (actinic)  fluxes  (horizontal  irradiances  arc  smaller  than  actinic  fluxes), 
and  their  effect  on  predicting  the  number  of  days  a  spherical  virus  can  survive  in  the  atmosphere. 
In  this  study,  we  compute  the  regression  model  only  for  MODTRAN  plane- irradiances.  Our 
regression  model  assumes  the  specific  shape  of  an  action  spectrum  that  is  responsible  for  the 
inactivation  of  the  viruses. 

MODTRAN  is  a  numerical  code  that  is  accepted  by  the  scientific  community,  and 
it  is  commonly  used  to  address  diffuse  sky  radiance,  multiple  scattering,  and  band-models  gas 
absorption.  MODTRAN  has  been  extensively  verified  and  validated  ( 4 ).  The  main  drawbacks  of 
MODTRAN  are  its  complexity  (the  user  must  be  an  expert  in  MODTRAN  to  be  able  to  use  it 
correctly)  and  its  application  is  time  consuming  (many  MODTR  AN  runs  arc  needed  to  compute 
U  V  daily  energy).  In  addition  (true  for  any  numerical  code),  the  output  of  each  run  from 
MODTR  AN  is  numbers;  thus,  intuitive  insight  for  the  effect  of  key  parameters  on  the  output  is 
lost.  To  overcome  MODTRAN  complexity,  wc  use  a  simple  nonlinear  regression  model  that 
gives  the  user  the  functional  dependence  of  UV  on  key  environmental  parameters  (altitude, 
albedo,  and  visibility). 

The  importance  of  a  regression  model  for  MODTRAN  results  cannot  be 
overstated.  The  size  of  MODTRAN  data  used  to  compile  the  regression  model  (Tables  1-5)  is 
~2  GB  (only  for  MODTRAN  Jlx  and  tp5  data  files),  and  it  takes  hours  of  central  processing  unit 
(CPU)  time  (on  a  fast  personal  computer).  We  acknowledge  that  with  greater  expertise  of 
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MODTRAN  execution,  data  storage  and  CPU  time  may  be  optimized.  We  also  acknowledge  that 
there  are  many  available  products  (5)*  that  ease  the  complexity  of  operating  MODTRAN  and 
may  reduce  CPU  time,  but  the  regression  model  proposed  here  is  virtually  instantaneous. 

The  novelty  and  usefulness  of  our  comprehensive  regression  model  is  that  it  ean 
“replace'”  MODTRAN  for  this  particular  problem  (inactivation  of  microbes  suspended  in  the 
atmosphere  by  solar  UV  radiation)  without  reducing  accuracy  significantly.  This  approach 
contrasts  previous  studies  where  MODTRAN  output  for  specific  limited  parameters  (e.g.,  UV 
transmission  [6],  water  temperature  retrieval  [7],  and  radiance  along  a  specific  flight  line  [<V]) 
were  fitted  with  a  regression  model.  These  types  of  regression  fits,  while  being  useful,  were  not 
substitutes  for  the  MODTRAN  vast  complex  output  that  our  regression  model  was  able  to 
capture.  Our  goals  are  to  demonstrate  that  a  comprehensive  regression  model  is  feasible  (for  the 
limited  UVB  wavelength  range  in  our  germicidal  application),  to  give  details  of  how  to  develop 
it,  and  to  encourage  others  to  attempt  this  course  of  action  for  other  MODTRAN  applications. 

This  report  is  organized  as  follows.  In  Section  2,  we  briefly  present  the 
terminology  and  geometry  of  the  radiative  transfer  for  virus  inactivation  in  the  presence  of  UV 
radiation  in  a  MODTRAN  atmosphere.  In  Section  3,  we  outline  a  simple-to-use  nonlinear  least- 
squares  regression  model  for  UV  radiation  fields,  sueh  as  total  daily  energy,  maximum  ambient 
flux  at  noon,  and  various  ratios  of  energy  to  flux,  that  are  key  components  for  estimating  virus 
inaetivation  in  theoretical  and  experimental  studies.  In  Section  3.1 ,  we  explore  the  range  of  the 
environmental  parameters  (altitude,  ground  albedo,  and  sea-level  visibility  distance)  for  the 
regression  model,  and  in  Section  3.2,  we  give  details  on  how  to  compute  the  regression  model 
coefficients  from  MODTRAN  simulated  atmospheres.  In  Section  4,  we  compare  the  results 
obtained  with  the  regression  model  with  MODTRAN  simulations  .  In  Section  4.1,  we 
demonstrate  that  the  radiation  fields  can  be  adequately  described  with  a  second  order  regression 
model.  Section  4. 1  also  provides  a  regression  model  for  several  atmospheres  for  a  35  latitude 
and  equinox  and  solstice  dates  in  the  form  of  relatively  simple  tables  and  equations  that  provide 
a  means  for  a  quick  estimate  of  UV  radiation  fields  and  inaetivation  time.  In  Section  4.2,  we 
attempt  the  sealing  of  various  radiation  fields  to  different  geographic  locations  (  latitude)  and 
dates  (day  of  the  year)  to  enable  the  use  of  the  regression  model  in  any  given  location  and  time 
of  the  year.  We  conclude  with  a  summary  in  Section  5. 


2.  RADIATIVE  TRANSFER  MODEL 

The  radiative  transfer  model  provides  a  measure  of  energy  conservation  for  the 
input  solar  UV  (at  the  top  of  the  atmosphere).  The  atmosphere  transforms  the  input  solar  photons 
into  one  of  the  following  types  of  photons: 

•  Scattered  photons  (reflection  of  photons  by  the  ground  is  technically 
considered  scattering) 

•  Absorption  photons 


*  Products  sueh  as  the  Electro-Optical  Systems  Atmospheric  Effects  Library  (EOSAEL)  and 
Ultraviolet  Transmission  and  Lidar  Simulation  Model  (UVTRAN). 
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•  Emission  photons  (emission  is  usually  negligible  in  UV  wavelengths) 

The  optical  properties  of  the  atmosphere  are  vertically  distributed  by  the  following  methods: 

•  Molecular  scattering  by  trace  gases  (Rayleigh  scattering  [9]) 

•  M1E  scattering  (9)  by  particulates  (aerosols) 

•  Absorption  by  molecules  and  aerosol  particles 

Scattering  is  not  restricted  to  one  interaction  (i.e.,  a  single-scattering  event) 
between  a  photon  and  a  single  molecule  (or  particle).  A  photon  may  be  rescattered  many  times 
(multiple  scattering,  namely  diffuse  light)  by  molecules  and  aerosols.  Aerosol  scattering  is  a 
directional  angular  scattering  where  the  angular  distribution  of  the  scattering  is  a  complicated 
function  of  the  chemical  composition,  size,  and  shape  of  the  aerosols.  The  molecular  absorption 
for  infinitesimal  monochromatic  radiation  is  described  by  the  Beer-Lambert  law  (9)  and  by  a 
band  model  (2)  for  radiation  with  a  given  spectral  width.  To  complicate  matters  further,  some  of 
the  scattered  photons  are  reflected  by  the  underlying  surface  (ground),  some  of  which  undergo 
further  scattering.  In  our  model,  we  assumed  spectrally  uniform  ground  reflectance  for 
convenience. 

We  divided  MODTRAN  parameters  into  spatial  and  temporal  parameters.  Spatial 
parameters  affect  the  composition  of  the  atmosphere  and  temporal  parameters  affect  the  time 
dependence  of  UV  radiation,  and  thus,  they  are  important  for  computing  the  daily  UV  energy 
incident  on  a  virus  (i.e.,  integration  of  flux  [irradiance]  over  time)  where  the  solar  UV  zenith 
angle  varies  as  a  function  of  time.  Spatial  parameters  are  determined  as 
y/  =  y/ (atmospheric  optical  properties)  ,and  a  temporal  parameter  is  set  at  0,  which  is  the  solar 
zenith  angle  as  the  function  of  the  day  of  year,  time  of  day,  and  geographical  latitude  and 
longitude  (in  this  work,  longitude  is  zero);  0=  0  is  the  zenith  angle  at  local  noon. 

The  notations  and  the  geometry  of  the  MODTRAN  atmosphere  were  summarized 
in  our  previous  paper  (/),  and  they  can  also  be  found  in  detail  in  the  MODTRAN  user  guide 
(10).  The  objective  of  the  present  study  was  to  compute  the  radiation  field  at  altitude  h  (above 
the  ground)  in  the  atmosphere.  The  spectral  flux,  F(Wlm2lnm),  per  wavelength  interval  v  at 
altitude  h  is  the  sum  of  the  attenuated  direct  flux  from  the  sun,  upward  diffuse  flux,  and 
downward  diffuse  flux.  All  fluxes  at  altitude  h  are  a  function  of  the  zenith  solar  angle  0  and  the 
optical  properties  of  the  atmosphere  y/.  We  stress  that  both  upward  and  downward  diffuse  fluxes 
(identified  in  the  equations  with  up  ]  and  down  ]  arrows)  are  affected  by  the  ground  albedo 
due  to  multiple  scattering  and  reflections  off  the  ground.  All  the  fluxes  (irradiances)  considered 
in  this  paper  are  plane  (horizontal)  irradiances  (i.e.,  irradiances  that  cross  a  horizontal  unit 
surface). 

A  virus  at  altitude  h  above  the  ground  is  illuminated  by  two  types  of  radiations: 

•  Direct  (unscattered)  solar  UV  flux  attenuated  by  the  atmosphere  along  the 
direction  0  to  the  virus  location 
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•  Diffuse  flux  due  to  the  complex  interaction  between  the  solar  UV 
radiation,  the  atmosphere,  and  the  ground 

The  total  flux  in  the  atmosphere  is  shown  by 


FT  (h,9,y/  ,v)  =  Fd,f  (h,9.<{/.v)  +  F*{h,0,y,v)  (1) 

where  superscript  T stands  for  total,  dif  for  diffuse,  and  dir  for  direct.  The  directional  upward  and 
downward  fluxes  at  altitude  h  were  combined  to  one  diffuse  flux: 

FJlf(h,o,i//,i’)  =  Fdlf(h,0,i//>v')+  F^r(h,0,i//,\’) .  The  portion  of  the  solar  UV  spectrum  efficient  in 

inactivating  viruses  lies  within  280-320  nm  and  it  is  generally  known  as  “action  spectrum” 

(//,  12).  We  computed  two  types  of  spectrally  integrated  fluxes  (eq  2):  Atmospheric  ambient 
flu \Fx(h,$,y/){W!m1)  and  Action  flux  Fax(h.0,<r)(lV/nr)  (this  flux  is  responsible  for  the  inactivation 
of  the  virus).  The  subscript  a  represents  action  and  the  superscript  x  is  reserved  for  total,  diffuse 
and  direct  radiation.  The  total,  diffuse  and  direct  spectrally  integrated  fluxes  are  shown  by 

320  nm 

F'{h,9,tf/)=  \F'(h,9,t//,v)dv 

l  (2) 

320nw  v  7 

F*(h,  9.i//)  =  jFx(h,O,<y,v)0Jv)d\’ 

„  280  nm 


where  is  an  action  spectrum  (/)*  that  gives  the  spectral  flux  fraction  (i.e.,  a  transmission 
function)  needed  for  virus  inactivation.  This  action  spectrum  is  believed  to  represent  the 
inactivation  sensitivity  of  diverse  viruses  belonging  to  viral  families  whose  members  cause 
severe  diseases  (//). 

The  action  spectrum  used  in  this  work  represents  viral  survival  at  various 
wavelengths  with  respect  to  the  survival  measured  after  exposure  to  254  nm  UV  radiation  (peak 
of  DNA  absorption  and  viral  damage).  Microbial  sensitivity  is  generally  expressed  by  the  D37 
survival  rate,  which  corresponds  to  the  irradiance  needed  to  produce,  on  average,  an  inactivating 
hit  per  virion,  which  corresponds  to  a  decrease  of  the  measured  survival  to  37%  of  the  viral 
amount  originally  irradiated.  Although  the  action  spectrum  was  built  using  D37  values  measured 
at  various  wavelengths,  other  survival  levels  (/)  can  also  be  employed  (e.g.,  the  survival  levels 
for  1  log'"  or  2  log1"  are  10%  and  1%,  respectively).  Our  regression  model  assumes  this  specific 
spectral  shape  of  the  action  spectrum;  a  different  spectral  shape  (e.g.,  for  a  different  virus)  will 
result  in  different  numerical  values  for  the  regression  model  coefficients  and  will  affect  the 
accuracy  of  the  regression  model.  The  virucidal  action  spectrum  was  determined  w  ith  data  from 
many  different  viruses,  where  the  target  is  the  nucleic  acid,  and  thus,  it  is  unlikely  that  a  virus 
will  show  a  different  action  spectrum  than  the  one  we  use. 

We  were  interested  in  estimating  the  daily  UV  energy  (a  full  day  exposure)  that  a 
virus  would  experience  at  a  specific  location,  even  though  the  virus  may  not  have  maintained  the 
same  location  from  sunrise  to  sunset  because  of  wind,  turbulence,  and  gravitational  settling.  For 


*  An  action  spectrum  is  illustrated  in  Figure  3  of  reference  (/). 
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computing  the  UV  inactivating  daily  energy,  we  needed  to  track  the  solar  angle  as  a  function  of 
time.  For  eaeh  point  in  time,  a  new  set  of  atmospheric  fluxes  was  computed  (eq  2).  The  solar 
zenith  angle#  =  ^(latitude,  day,  hour)  (13)  as  a  function  of  latitude,  day  of  year,  and  time  of  day, 
is  given  by 


cos(0)  =  sin(/a//7ude)sin(J)  +  cos{latitude)cos(S)cos,(ha) 

d  =  0.00691 8  -0.39991 2 cos(y)  +  0.070257 sin(^) -  0.006758  cos(2/)  +  0.000907 sin(2^)-... 

0.002697  cos(3^)  +  0.001 48  sin(3^) 

2^((/m'  -  1  +  (/tow  -  12)/  24 

7  ~  365 

,  31 
ha  = - x 

180 

60 hour  +  229.18(75  x  10  *  +  0.001868cos(^)  -  0.03 2077 sin(^ )-  0.0 146 15 cos(2/)  -  0.040849 sin(2/)) 

4 

,  ..  ,  2  x  sunrise  180.  .  v 

daylight  - -  ;  sunrise  = - |arccos(-  lan (latitude)  lan(c>))[ 


(3) 


80) 


where 

<5  is  the  solar  declination  angle  (radians) 

ha  is  the  solar  hour-angle  (radians) 

y  is  the  fractional  year  in  radians 

day  is  the  day  of  the  year  (1  January  is  day  1) 

hour  is  the  hour  in  the  day  (0  to  24) 

sunrise  is  the  hour  angle  (degrees)  for  sunrise 

daylight  is  the  daylight  time  (hours)  between  sunrise  and  sunset 

In  cq  3,  time  of  day  is  local  solar  time  (hour  =  12  noon;  the  sun  is  at  the  zenith 
and  hour  can  be  given  as  a  fraction  [e.g.,  hour=  1 1.5]).  The  sunrise  and  sunset  times  are  the 
times  for  which  the  solar  zenith  angle  is  at  the  horizon  (i.e.,#  =  90° ).  The  sunrise  and  sunset 
times  are  symmetric  with  respeet  to  noon.  Therefore,  the  length  of  daylight  is  easily  computed  as 
(sunset  -  sunrise)!  15  =  2  x  sunrise!  15 ,  where  the  constant  1 5  is  the  value  converted  from  degrees  to 
hours  (the  earth  rotates  at  an  angular  speed  of  1 5°/h).  We  ean  write  the  time-dependent  solar 
angle<9as0(/7,/),  where denotes  the  geographical  latitude  and  the  day  of  the  year, 

»  t)  =  (latitude,  day  of  year) ,  and  t  is  the  time  (in  hours)  within  the  day.  With  these  notations,  the  daily 

energy  E(J  lm2)  at  altitude  h  and  locatioa^day  7  is  the  integration  of  fluxes  (eq  2)  over  time  of 
«  day  (/).  This  equation  is  given  by 
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sunset 

Ex(hjj,y/)  =  J  Fx(h,0(rjJ),i//)Ji 

sunrise 

sunset 

<  E'a{fut],\}/)  =  \F*{h,6(t],t),yt)dt  > 

sunrise 

rj  =  q  {latitude^  day  of  year) 


(4) 


where 

Ex(h,q.ty)  is  the  atmospheric  daily  energy 

E'a(h,t],y/)  is  daily  action-energy 

The  portion  of  the  diffuse  daily  action-energy  relative  to  the  total  is  variable 
(from  85%  near  the  ground  to  53%  at  5  km  for  U.S.  1976  standard  atmosphere  in  spring-summer 
with  rural  aerosols  on  21  March  and  latitude  of35°N)  (/).*  It  is  a  function  of  visibility,  albedo, 
and  altitude  (increasing  with  albedo,  decreasing  with  visibility  and  altitude).  Many  monitoring 
institutions  measure  and  report  the  ambient  flux  only  at  local-noon  when  the  sun  is  at  the  highest 
point  in  the  sky.  These  local-noon  measurements  correspond  to  the  maximum  magnitude  of 
spectrally  integrated  ambient  flux  where  the  ambient  flux  attains  its  maximum  magnitude: 

F'  ->  .  The  maximum  value  of  ambient  tlux  is  given  by 

320 wm 

F^x(h,t},y/)  =  F'  (h,0(i},t  =  noon),y/)  =  \  F1  (hy0{i]j  =  noon\y/  ,v)dv  (5) 

280rtm 


3.  METHOD  OF  SOLUTION 

For  predictions  of  virus  inactivation,  the  theoretician  who  models  the  radiative 
transfer  interaction  and  the  experimentalist  who  measures  the  UV  flux  (or  daily  energy)  may  be 
interested  in  different  quantities.  The  theoretician  may  be  interested  in  the  daily  action 
energies,  E*  ,  Ef  (h,rj,t//) ,  and  Efr(h,i]yy/) ,  as  functions  of  atmospheric  parameters  ^and 

location  (h.rj) .  The  experimentalist  who  measures  the  ambient  atmospheric  fluxes  at  noon  may 
want  to  seale  the  noon-time  measurements,  F^x(h, /;,(/),  F^(h,ij,iy) ,  and  F^(/»,;;,^) ,  to  daily 
action  energies. 

Computing  daily  action  energies  with  MODTRAN  is  complex  and  time- 
consuming.  For  each  time  increment  t  and  atmospheric  condition  y  in  cq  4  (e.g.,  surface  albedo 
or  meteorological  visibility  range),  a  new  MODTRAN  run  must  be  made  and  numerical 
integration  of  the  results  must  be  performed.  Our  objective  in  this  study  is  to  develop  a  relatively 
simple  regression  model  that  reproduees  MODTRAN  results  and  eqs  2,  4,  and  5  within 
acceptable  accuracy  without  intensive  computations.  To  the  best  of  our  knowledge,  a  regression 
model  has  not  been  attempted  successfully  on  MODTRAN.  However,  in  this  study  (it  may  be 
due  to  the  spectral  shape  of  the  action  speetrum  and  the  relatively  narrow  wavelength  range),  we 


*  Illustrated  by  Figure  5  in  reference  (/). 
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were  pleasantly  surprised  to  see  that  reasonable  accuracy  in  the  regression  model  was 
accomplished  (Section  4). 

Our  simple  regression  model  addresses  nine  quantities:  three  daily  action  energies 
(total,  diffuse,  and  direct)  and  six  ratios  given  by 


— ■  (s)  =  total  daily  action- energy  normalized  to  total  ambient  flux  at  noon 

JJh’hY) 

piiif  f  \ 

— J- — 'JM—  (s)  =  diffuse  daily  action- energy  nonnalized  to  diffuse  ambient  flux  at  noon 

FjLfa’h'P) 

(s)  =  total  daily  action- energy  normalized  to  direct  ambient  flux  at  noon 

— 1 7 — =  fraction  of  diffuse  daily  action- energy  in  total  daily  action- energy 
Ea(h,rj,y/) 

!LA!hIb¥l  =  ratio  of  direct  daily  ambient  energy  (w/o  effect  of  virus  action  spectrum) 
to  totaldaily  action- energy 

— ~ — — --  =  fraction  of  ambient  diffuse  flux  at  noon  in  total  ambient  flux  at  noon 

ET2(h,r},y/)  =  total  daily  action- energy(J/nr) 

Ed/ {h,T],y/)  =  diffuse  daily  action -energy(J/m) 

E‘‘ir(h,T},y/)  =  direct  daily  action- energy(J/m2) 


(6) 


All  of  the  quantities  in  eq  6,  except  £ar  and  Ef ,  are  new  and  were  not  included  in 
our  previous  paper  (1).  For  practical  reasons  (to  be  explained  in  Section  3.2),  we  chose  to 

compute  the  inverse  (reciprocal)  of  and  (see  Section  4,  Tables  1—4). 

Ea(h,rj,</) 

Several  other  useful  quantities  can  be  derived  from  the  nine  quantities  in  eq  6.  For  example,  the 
ambient  atmospheric  fluxes  (i.e.,  total,  direct,  and  diffuse  ambient  fluxes  that  illuminate  the  virus 
at  altitude  /;,  without  the  effect  of  the  virus  action  spectrum)  and  the  ambient  atmospheric  direct 
(unscattered)  daily  energy  can  be  calculated  by 


Fl{hw)  =  E?{Kn,¥\ 


'Jif  t 


FlV>W) 


Fi FZ(h,r),y) 

rdir,,  X  rT / 1  , 

Ea(h,n,y) 


(7) 
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With  our  regression  model,  we  can  retrieve  13  quantities  of  potential  biological 
utility  as  shown  by  eqs  6  and  7.  We  caution  the  reader  that  with  the  regression  model  coefficients 
(given  in  Tables  1~4  for  eq  6),  cq  7  results  show  an  amplification  of  errors  (e.g.,  the  relative 
error  in  a  calculated  ratio  x/y  where  .v  and  y  suffer  1 0%  error  is  20%  for  uncorrclatcd  errors), 
which  lowers  the  accuracy  obtained  by  eq  6. 

After  carefully  considering  the  contribution  of  various  parameters  to  virus  » 

irradiation  in  the  atmosphere,  three  parameters  were  selected  for  the  regression  analysis: 

•  The  altitude  of  viruses  above  the  ground  h 

•  The  surface  albedo  (reflectivity)  of  the  land  underneath  the  airborne 
virsues 

•  The  visibility  range  (inversely  proportional  to  aerosol  number  density) 

The  relative  humidity  (RH)  is  also  highly  variable,  but  its  effect  is  accounted  for 
in  MODTRAN  simulations  by  the  visibility  distance  parameter.  We  pursued  a  regression  model 
in  the  form  \ih,atb.vis\q,y)  that  captured  the  parameter  space  (altitude,  albedo,  and  visibility)  for 
each  of  the  nine  quantities  (six  ratios  and  three  daily  action-energies)  at  a  given  geographic 
location  and  day  of  the  year  (q)  and  atmosphere^). 

The  effect  of  the  vertical  transport  of  viruses  due  to  air  currents  could  be 
incorporated  in  the  model  (in  principle)  by  computing  a  regression  model  as  a  function  of 
altitude,  visibility,  and  albedo  for  each  hour  and  summing  the  altitude-dependent  energy  within 
time  span  t /  to  t;,  if  the  whereabouts  (altitude)  of  the  virus  is  known.  We  chose  only  to  compute 
the  energy  within  daylight  (sunrise  to  sunset)  to  maintain  a  relatively  simple  approach.  We  note 
that  the  falling  velocity  of  a  small  spherical  particle  in  the  atmosphere  is  very  small;  a  particle 
<1  pm  diameter  (the  size  of  an  agglomerate  of  viruses  is  on  the  order  of  microns,  the  size  of  a 
single  virus  is  in  the  range  of  0.01  to  0.1  pm)  with  a  density  of  1  g/m  will  fall  in  calm  air  (no  air 
turbulence  and  no  vertical  air  currents)  at  a  terminal  gravitational  settling  velocity  of -3x10  5 
m/s.  Thus,  within  24  h,  a  particle  will  fall  less  than  3  m. 

3.1  Parameter  Space  for  Regression  Model 

Altitude. 

# 

We  restricted  the  altitude  h  above  the  ground  to  a  range  of  values  0<  h  :  ?  km .  These 
values  encompass  most  scenarios  of  interest  for  the  potential  location  of  viruses  in  the 
atmosphere.  In  the  numerical  grid  for  /;,  we  used  more  values  near  the  ground  in  order  to  capture 
the  complexity  of  the  radiative  transfer  model  in  the  boundary  layer  (where  most  of  the  ambient 
atmospheric  aerosols  reside).  Twenty  values  were  chosen  for  h:  0,  0.005,  0.01 , 0.02,  0.05,  0.075, 

0. 1 , 0. 1 5,  0.2,  0.3,  0.4,  0.5,  0.75,  1 .5,  2,  2.5,  3,  4,  and  5  km. 
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Sea-level  visibility  ranee. 


r 


The  sea-level  (surface)  visibility  meteorological  range,  vis,  in  MODTRAN  is 
based  on  1926  Koschmieder’s  theory  (14),  where  visibility  is  given  as  a  function  of  the  surface 
aerosol  extinction  at  wavelength  550  nm: 


vis  (km)  = 


_ T912 _ 

aerosol  extinction  (km  ')+ 0.01 159  (Am1) 


(8) 


where  0.01 159  is  the  molecular  scattering  coefficient  at  550  nm  and  3.192  is  log(2%)  for  the 
assumed  2%  brightness  contrast  with  which  the  visibility  range  is  defined. 

Equation  10  shows  that  the  aerosol  loading  in  the  atmospheric  boundary  layer 
(i.e.,  surface  aerosol  extinction  in  the  denominator)  determines  the  visibility  distance.  In  our 
study,  we  used  this  (MODTRAN)  definition  for  visibility  range.  Later  studies  (15)  revised  the 
relationship  between  reported  visibility  ranges  (by  observers  in  meteorological  stations)  and  an 
aerosol  extinction  coefficient  to  be  obsen-ed  vis(km)  =  (1 .9/3.9 1 2)  *  vis  =  vis  12  because  visual  targets 
used  by  the  observer  to  assess  visibility  distance  are 

•  Not  black 

•  Frequently  too  small  in  angular  size 

•  Located  only  at  quantized  distances 

This  modified  visibility  range  (15)  was  used  (16)  to  estimate  the  distribution  of 
continental  surface  aerosol  extinction.  It  was  found  that  the  extinction  coefficient  in  the  eastern 
United  States  was  relatively  modest  (0. 1-0.2  km'1).  The  observed  visibility  range  was 
9.5-19  km,  which  corresponds  to  a  MODTRAN  visibility  of  19-38  km. 

In  UV  wavelengths,  a  RH  >  ~75%  alters  the  optical  properties  of  aerosols  (such 
as  the  single-scattering  albedo)  and  thus  modifies  aerosol  extinction.  Most  of  the  ambient 
aerosols  reside  in  the  atmospheric  boundary  layer  (0-2  km),  where  the  effect  of  RH  is 
pronounced.  When  the  visibility  distance  is  set  in  MODTRAN  (cq  8)  the  number  density  of  the 
aerosols  (i.e.,  aerosol  loading)  is  sealed  to  produce  the  appropriate  extinction.  Thus,  when  the 
RH  in  a  MODTRAN  atmosphere  is  modified  (as  an  input),  the  effect  of  RH  on  the  aerosols  is 
mostly  cancelled  by  rescaling  the  extinction  (that  was  altered  by  RH)  with  the  visibility 
parameter.  For  this  reason  (verified  in  our  simulations),  we  chose  not  to  have  RH  as  a  parameter 
in  the  regression  model.  In  the  regression  model,  wc  restricted  MODTRAN  visibility  (cq  8)  to 
the  range  1 5  km  s  vis  s  too  km ,  which  spans  the  observed  visibility  range  of  7.6-69  km.  Eleven 
values  were  chosen  to  span  the  visibility  space:  15,  20,  25,  30,  40,  50,  60,  70,  80,  90,  and  100 
km. 
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Surface  albedo. 


The  measured  UV  albedo  values  of  most  natural  surfaces  are  smaller  than  10%. 

The  albedo  (MODTRAN  surref  parameter)  was  assumed  to  be  spectrally  uniform  between  280 
and  320  nm.  Several  spectral  albedo  values  (MODTRAN  internal  database  in  [2])  at  200  nrn 
(and  values  in  parentheses  at  500  nm)  are 

•  Grassland:  0.025  (0.065)  » 

•  Mixed  forest  (consisting  of  70%  broadleaf  and  30%  pine):  0.025  (0.036) 

•  Urban  (commercial):  0.026  (0.036) 

•  Deciduous  needle  forest  (broadleaf-pine):  0.017  (0.032) 

•  Closed  shrubs  (pine  brush):  0.012  (0.049) 

•  Tundra:  0.032  (0.07) 

•  Savanna  (grass-scrub):  0.0 1 9  (0.066) 

•  Woody  savanna  (soil-grass-scrub):  0.021  (0.06) 

•  Crop  mosaic  (mixed  vegetation):  0.007  (0.017) 

For  materials  tested, 

•  The  albedo  levels  of  old  grass,  dead  grass,  burnt  grass,  and  maple  leaf  at 
300  nm  were  0.07,0.1 , 0.02,  and  0.02,  respectively. 

•  The  albedo  levels  for  desert  sand  (mixture  of  roeks  and  silt  sand)  at  275, 

300,  and  325  nm  were  0.12,  0.13,  and  0.14,  respectively. 

•  The  albedo  level  for  ocean  water  at  300  nm  was  0.027. 

In  the  UVB  region  (v  <  3 1 5  nm).  the  measured  albedo  (77)  was  as  small  as 
0.016-0.01 7  over  vegetation,  0.04-0.05  over  bare  fertile  soil,  and  0.07-0. 10  over  concrete 
(autobahn,  Germany).  The  albedo  over  dry  bright  sand  (e.g.,  ocean  beach)  was  0.14  at  300  nm. 

3.2  Second  Order  Regression  Model 

In  this  section,  we  give  details  of  how  to  compute  the  regression  model  from 
MODTRAN  results  to  enable  the  reader  to  construct  regression  models  for  other  atmospheric 
conditions  and  radiative  transfer  scenarios.  We  experimented  with  different  options  for  a  « 

regression  model,  such  as  a  different  polynomial  order  with  each  variable  (h,alb,vis) ,  to  optimize 
the  regression  model  fit  to  the  simulated  data  (MODTRAN  output)  for  a  given  atmospheric  A 

model  </7,y/) .  We  searched  for  the  lowest  order  polynomial  that  showed  reasonable  results  and 
ensured  numerical  stability  of  the  regression  model.  In  the  regression  model,  we  treated  each 
variable  independently  and  did  not  allow  interaction  (cross-terms)  between  variables  (e.g., 
h x  alb  or  vis  x  alb )  to  preserve  physical  meaning  as  much  as  possible.  We  finally  chose  a  second 
order  polynomial  y(h,aib,vis)  with  seven  coefficients,  h’i  to  w-j  (determined  with  least-squares  fit  to 
MODTRAN  data),  for  given  atmospheric  conditions  (w)  by 
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y(h,alb,  vis  \  T],y/)  =  w,  +  h 2h  +  v.\h:  +  w<a!b  +  wsa/b2  +  w6vis  +  w-,vis2 


(9) 


•r 


where  h  is  altitude  in  kilometers,  visibility  is  in  kilometers,  and  the  albedo  is  a  fraction  (e.g., 

1 0%  reflectivity  is  an  albedo  of  0, 1 ).  We  used  eq  9  for  each  of  the  nine  quantities  listed  in  eq  6. 
However  (for  reasons  of  numerical  stability,  discussed  below),  we  chose  to  do  a  regression 
model  y(h,aib,vis)  for  the  inverse  of  E\  /  F*'x  and  Ef  /  £ar .  These  two  quantities  are  easily  retrieved 

by  - .  - ! - and  —  - 1 - -  The  other  seven  quantities  (four  ratios 

F£[h,n,v)  yfh,alb,vis)  ETa(h,ih¥)  y(h,alb,vis) 

and  three  daily  action-energies)  are  shown  directly  with  the  regression  models 

=y(h,alb,vis\rj,y/)  and  £ar(A,//,^)  =  y(h,alb,vis\i],y/) .  Assuming  that  the  energy  needed  to 

inactivate  virus  X  to  a  given  survival  level  is  E,  (J/nE)  (/,  11),  using  our  regression  model,  we 
can  estimate  the  maximum  number  of  radiation-days,  ,  needed  for  inactivating  virus  X  at 
altitude  h  in  a  given  atmospheric  model,  visibility  condition,  and  ground  albedo  by 

Hj,n,(b,alb,vis)  = - — - = - - - — - ; - 7 

y(h,alb,vis)  wl+w2h+wih  +  wAath  +  w,alb  +wivis  +  w1vis 

(10) 


where  yih,aib,vis)  refers  to  any  of  the  last  three  quantities  in  eq  8.  Equation  10  provides  the 
maximum  number  of  radiation-days  because  the  plane  irradiances  underestimate  the  true 
spherical  irradiances  (see  [/]  for  a  discussion  on  the  difference  between  plane  and  spherical 
fluxes).  The  complexity  of  the  effects  of  altitude,  albedo,  and  visibility  on  inactivation  time  was 
previously  demonstrated  and  discussed  (/)*,  where  the  number  of  radiation-days  needed  for 
inactivation  is  1  to  2  radiation-days  in  winter  (at  latitude  35°  N,  on  21  December). 

To  solve  eq  10  for  the  seven  coefficients,  we  constructed  the  matrix  equation 
Y  =  A  w  ,  where  Y  is  a  vector  constructed  from  MODTRAN  output 

y(h,alb,vis\n,y/)={MODTRAN results)  for  all  (1760)  combinations  (h, alb, vis) ,  A  is  a  matrix,  and  w  is  a 
vector  of  the  unknown  seven  regression  coefficients.  The  matrix  A  and  vectors  Y  and  w  are 
shown  by  eq  11. 


*  Illustrated  by  Figure  6  in  reference  (/). 
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To  solve  for  h\  we  performed  inverse  operation 
>v  =  A  ^ 


(12) 


where  the  nonsquare  matrix  A  J  is  computed  by  singular  value  decomposition  (/<V)  using  seven 
eigenvectors  and  eigenvalues  (because  the  rank  of  the  matrix  A  is  rank(/l)  =  7)  or  by  a  Gaussian 
elimination  algorithm  with  partial  pivoting.  The  error  magnification  (i.e.,  the  stability  of  the 
regression  model  to  input  errors  in  Y  in  the  solution  w)  is  proportional  to  the  condition  number 
(ratio  of  the  largest  to  the  smallest  eigenvalues  of  A).  Due  to  the  large  dynamic  range  in  A  (/; 
ranges  from  0  to  5  km,  albedo  ranges  from  0  to  0.15,  and  visibility  ranges  from  15  to  100  km),  to 


achieve  numerical  stability,  we  normalized  the  variables,  h 


alb 


albm 


,  and  vis 


v/.s- 


in  the  matrix  A.  As  a  result  of  the  normalization,  we  had  to  rescale  the  solution  "  *  A  '  so  that 
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remained  unchanged)  to  obtain  a  stable  regression  model  (eq  9).  With  this  normalization 
procedure,  the  condition  number  improved  (i.e.,  it  decreased)  from  236,410  to  33.3.  The 
maximum  error  magnification  equals  5.7  =  >/33  and  the  mean  relative-error  magnification  (19) 


was 


|(m] /jw|  ^/max  eigenvalue/ min  eigenvalue 


,  where  \\  denotes  a  vector  norm,  <*>•/*  is  the  relative 


error  in  tv,  and  d\/Y  is  the  relative  error  in  MODTRAN  output.  The  mean  error  magnification  for 
the  seven-element  veetor  tv  is 0.82  =  1/33/7.  This  value  is  small  enough  to  ensure  stability  and 
accuracy  of  the  solution  tv  for  the  regression  model  given  in  eq  9. 


4.  RESULTS 

We  computed  a  regression  model  (Tables  1-4)  for  MODTRAN  atmospheres 
(MODTRAN  parameter  model  =  2,  3,  6  ),  which  showed  1976  U.S.  standard  and  mid-latitude 
atmospheres  for  two  types  of  ambient  aerosols.  These  atmospheres  represented  rural  and  urban 
environments  (MODTRAN  parameter  ihaze  was  either  1  or  4;  please  note  that  MODTRAN 
parameter  vis  alters  the  concentration  [number  density]  of  aerosols  in  the  atmospheric  boundary 
layer,  as  set  by  ihaze,  to  eomply  with  the  desired  visibility  distance).  For  each  atmosphere,  we 
evaluated  the  performance  of  the  regression  model  during  two  seasons:  spring-summer  and  fall- 
winter  (MODTRAN  parameter  iseasn  was  either  1  or  2),  and  four  dates  (equinoxes  and  solstices, 
21  Mareh,  21  June,  21  September,  and  21  December).  The  concentration  and  altitude  profiles  of 
all  the  gases  were  set  internally  by  MODTRAN  default  parameters.  In  the  computations,  wc 
assumed  a  elear  sky  and  a  speetrally  flat  albedo.  Beeause  a  large  portion  of  the  world  population 
lives  within  30  to  40°  latitudes  (20),  we  ehose  a  latitude  of  35°  for  our  calculations.  In  section 
4.2,  wc  show  our  scaling  (Table  5)  of  the  regression  model  to  other  latitudes  (25,  35,  and  45°) 
and  times  of  year  (1  February,  1  May,  1  August,  and  1  November). 

4. 1  Regression  Model 

Rational  for  a  seeond  order  regression  model. 

The  Figure  supports  the  selection  of  the  second  order  regression  model  described 
by  eq  9  with  altitude,  albedo,  and  visibility  as  variables.  The  attenuation  in  the  atmosphere  is  an 
*•  exponential  proeess  with  optical  depth  (e.g.,  Beer’s  law  for  the  direct  component  of  the  solar 

UV)  that  is  nonlinear  with  altitude,  as  is  shown  in  the  curvature  (steeper  at  lower  altitudes)  of  the 
direct  daily  action-energy  (Figure  c).  The  nonlinear  effect  of  altitude  is  also  evident  in  the 
eurvature  of  the  surfaces  in  Figure  b  and  in  the  nonuniform  spacing  of  the  surfaces  in  Figure  a 
(e.g.,  the  spaeing  between  the  1  and  3  km  surfaces  is  larger  than  the  spaeing  between  3  and 
5  km).  This  observation  led  us  to  choose  a  regression  model  with  second  order  polynomial  with 
altitude,  y(h,alb,vis\rj,y)<x  w^+  Wyh1  in  eq  9.  The  nonlinear  effect  of  visibility  is  evident  in  the 
curvature  of  the  surfaces  in  Figure  a  for  the  diffuse  daily  action-energy,  in  Figure  c  for  the  direct 
daily  action-energy,  and  in  Figure  d  for  the  total  (diffuse  +  direct)  daily  action-energy  where  the 
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polarity  of  the  curvature  is  concave  near  the  ground  and  convex  at  higher  altitude.  As  altitude 
increases,  the  curvature  of  the  surfaces  decreases. 

The  nonlinearity  with  respect  to  visibility  led  us  to  choose  a  regression  model 
with  the  second  order  polynomial  with  visibility  y(h,alb,vis\r),y/)  x  w6vis  +  w7»-/.vJ .  The  almost  linear 
effect  of  albedo  (linear,  but  with  a  nonzero  offset)  is  shown  in  Figure  a  for  the  diffuse  and  in 
Figure  d  for  the  total.  Even  though  the  effect  of  albedo  is  nearly  linear,  we  chose  a  second  order 
polynomial  with  albedo  \ih,alb,\is\rj,y/)  x  wAalb  +  w<alb2  for  two  reasons: 

(1)  The  ground-refleeted  direct  flux  is  proportional  to  albedo,  hence 

y(h, alb ,  vis  \  rf%\y)cc  w4alb  . 

(2)  T  he  contribution  of  diffuse  flux  due  to  multiple  reflections  off  the  ground 
is  proportional  to  higher  order  effects  in  albedo. 

At  short  wavelengths  (e.g.,  a  visible  wavelength  region  where  thermal  emission  is 
negligible),  the  diffuse  flux  is  proportional  to  aIb/(\-albxS)(2]),  where  S  is  the  atmospheric 
spherical  albedo  (the  portion  of  radiation  that  is  reflected  from  the  atmosphere  back  to  the 
ground  as  a  result  of  isotropic  illumination  from  the  ground).  A  second  order  Taylor  series 
expansion  of  alb  /(I  -  alb*  S )  is  shown  by  alb/(  1  -  alb*  S)  =  alb  +  S  x  alb 2 .  Thus,  the  second  order 
polynomial  with  albedo,  y(h,alb'Vis\n<v)'*  wsalb2 .  is  due  to  the  linear  effect  of  the  ground-reflected 
direct  solar  photons  and  the  quadratic  effect  for  diffuse  photons. 


20 


Figure.  Rationale  for  using  a  second  order  regression  model  for  altitude,  visibility,  and  albedo 
demonstrated  with  1976  U.S.  standard  atmospheric  model,  at  spring-summer  with  rural  aerosols, 
at  latitude  35°  on  21  March,  (a)  Daily  diffuse  action-energy  Eff(h, albedo, visibility  )  as  a  function 
of  visibility  range  and  ground  albedo  for  four  altitudes  above  ground  (in  ascending  order);  0 
(bottom  surface),  1, 3,  and  5  km  (top  surface),  (b)  Daily  diffuse  action-  energy 
E*f  (h. albedo,  visibility)  as  a  function  of  altitude  and  visibility  for  three  values  of  ground  albedos; 
0.12  (top  surface),  0.06  (middle  surface),  0  (bottom  surface),  (c)  Direet  action-  energy 
E^(h, visibility)  as  a  function  of  altitude  and  visibility  (the  direct  component  is  independent  of 
albedo),  (d)  Daily  total  (diffuse  +  direct)  action-energy  E^h, albedo, visibility) as  a  function  of 
visibility  range  and  ground  albedo  for  four  altitudes  above  ground  (in  ascending  order);  0 
(bottom  surface),  1 , 3,  and  5  km  (top  surface). 

Regression  model  for  various  atmospheres. 

Tables  1—4  show  the  application  of  the  regression  model  for  1976  U.S.  standard 
atmosphere  and  for  mid-latitude  atmosphere,  with  rural  and  urban  aerosol  models  at  the 
four  days  of  equinox  and  solstice  at  latitude  35°.  For  21  March  and  21  June,  we  used  a  spring- 
*•  summer  seasonal  model,  and  for  21  September  and  21  December,  we  used  a  fall-winter  seasonal 

model.  The  season  (parameter  iseasn  =  1  and  2  in  MODTRAN)  only  affected  the  aerosol  profile 
in  the  troposphere  (2-10  km)  and  stratosphere  (10-30  km).  There  are  two  mid-latitude 
atmospheric  models  (parameter  model  =  2  and  3  in  MODTRAN);  mid-latitudc-summcr  and 
mid-latitude-winter.  For  21  September  and  21  December,  we  used  (Tables  3  and  4)  the 
atmospheric  model  mid-latitude-winter,  and  for  21  March  and  21  June,  we  used  the  atmospheric 
model  mid-latitude-summer.  The  atmospheric  model  set  the  altitude  distribution  of  temperature. 
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pressure,  and  concentration  profile  of  trace  gases  in  the  atmosphere  (7).*  The  effect  of  the 
atmospheric  model  and  season  can  be  assessed  by  employing  the  regression  model  for  different 
choices  of  atmospheric  models  and  seasons  (7).  The  regression  model  coefficients  were 
computed  for  all  nine  quantities  (eq  6)  over  all  parameter  spaces  (h, alb, vis) .  The  mean  absolute 
error  and  the  maximum  absolute  error  are  given  in  the  tables  for  each  of  the  nine  quantities.  The 
mean  (average)  of  the  absolute  error  is  the  absolute  percentage  difference  between  the  regression 
model  in  eq  9  and  MODTRAN  output  (regarded  as  “truth”)  averaged  over  all  the  parameter 
spaces:  all  altitudes  (0<h<5km),  albedos  (0<a//><0.15),  and  MODTRAN  visibilities 
( 15  <  vis  £ lOOAw).  The  maximum  absolute  error  computed  for  a  probability  p  1 0.9  gives  the  error 
for  which  90%  of  all  y(A,  a/6,  v&)  regression-model  values  fall  within  MODTRAN  truth  values: 

\trmr\ 

max(jt77wj)  =  j  Probabiliiy([  error  |)  c7(|  error  I)  =  0.9 .  Closer  inspection  ol  the  coefficients  in  the  tables 

0 

reveals  that  because  e\  =Eff  +  ej"  ,  the  regression  coefficients  for  Efr  are  (within  round-off 
errors)  the  differenee  between  the  regression  coefficients  for/rj  and£^ .  The  tables  demonstrate 
(see  second  to  last  columns  in  each  of  Tables  1-4)  that  the  regression  model  shows  a  very  close 
approximation  to  MODRAN  results.  The  worst  accuracy  is  ~6%  for  radiation  fields  Ef  and 
F^/El  at  Dec  21.  For  all  other  quantities,  the  error  is  <5%,  and  in  most  cases,  it  is  <2%.  This 
agreement  appears  excellent,  especially  in  contrast  to  errors  in  microbiological  laboratory 
experiments  where  the  survival  level  of  viruses  to  UV  radiation  is  usually  measured  within  an 
accuracy  no  better  than  ±o.5iog,0  units  ( 1  log10  is  a  survival  level  of  10%).  We  remind  the  reader 

that  we  chose  to  do  a  regression  model y{h,aib,vis)  for  the  inverse  (reciprocal)  of  ETa  / Fn*[  and 
£‘0*  I FT.  •  Therefore,  Ef  / ETa  is  shown  (Tabic  1,21  March)  as 


_ , _ - _ , _ , — rand  the 

El  (h,  alb,  vis)  1.21  +  0.0647//  +  0.00387//'  -  0.75 1  alb  +  O.HSalh2  +  0.00696  vis  -  3.81  x  10  vis 

quantity  (eq  7)  F^(h,alb,vis)  is  easily  obtained  from  (Table  I,  21  March)  by 


F^(.h,alb,vis)  = 


55.7  +  17,6//  -0.923 h2  +  71 .2 alb  +  28.7  allE  +  0.  t48  vis  -0.000888 vis' \ 

26. t  +  0.979//  -  0.052//2  - 1.87  alb  +  0.066a/b2  +  0.0 17  vis  -  7.08  x  to  5 vis2 


where  h  and  vis  are 


in  kilometers  and  the  albedo  is  a  fraction  (0  to  1). 


4.2  Location  and  Time  Scaling  of  Regression  Model 

The  regression  model  in  Tables  1-4  was  computed  for  a  specific  latitude  (35°) 
and  day  of  the  year  (equinox  and  solstice  days).  Scaling  the  model  to  other  locations  (latitudes) 
and  times  (dates)  of  interest  should  save  enormous  effort  in  computing  new  regression 
coefficients.  Scaling  the  results  is  very  challenging  because  the  solar  zenith 
angle,  0(t, latitude, day  of  year)  ,(eq  3)  varies  with  location  and  time  of  day.  Hence,  the  input  UV 
solar  flux  at  the  top  of  the  atmosphere  (the  engine  for  all  subsequent  scattering)  is  proportional  to 
cos(0)  and  for  each  time  of  day,  we  have  a  radiative  transfer  model  with  a  new  input  solar  flux. 
The  zenith  angle  also  determines  slant-path  attenuations  within  the  atmosphere,  and  the 
interaction  with  the  altitude-dependent  gases  and  aerosols  is  vastly  different  as  a  function  ofo. 


i 


*  See  Figure  2  in  reference  (7)  for  the  differences  between  the  1976  U.S.  standard  and 
mid-latitude  summer  and  winter  atmospheres. 
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The  daily  energy  increases  with  the  length  of  daylight,  but  because  the  length  of  daylight  is 
linked  to  the  solar  angle,  the  relationship  between  daily-energy  and  length  of  daylight  is  far  from 
simple.  Nevertheless,  we  attempted  location/time  scaling  and  were  able  to  produce  reasonable 

Et  Edif 

results  for  six  quantities:  — f-;— '—■.-sfL\ET0\E?',E?  .  We  found  a  simple  scaling  with  latitude/day 

h  b  b 

mix  max  max 

as  a  function  of  solar  zenith  angle  and  length  of  daylight.  The  scaling  (for  any  given  atmospheric 
model)  from  a  reference  latitude  and  date  to  latitude  and  date  is  given  in  eq  16.  The  reference 
a  latitude  and  date  were  chosen  to  be  35°  and  2 1  March  for  the  given  atmosphere.  Although  we 

could  not  establish  a  physical  explanation,  we  noted  that  powers  of  cosine  angle  (e.g.,  second, 

,  third,  and  fourth  powers,  used  in  our  scaling)  are  common  in  determining  irradiance  produced  by 

a  diffuse  source  (22) 


y(h,alb ,  vis,  latitude,  day) 
y(h,  alb ,  vis , lat  -  35° ,  day  -  Mar  2 1) 


cos' [0(1  at,  day)]  [  E!(h,alb,  vis)  Ed,f  (h,  a  lb,  vis) 

- - - - -  for  v  =  <  — - -  and  - 

cos2  [0(35”,  Afor  21)]  [F^(h,alb,vis)  FZ(*>,alb,vis) , 

cos[0(35',  Afor  2 1)]  x  day!ight(lat,da\)  F*l(h,alb,vis) 

cos[0(lat,day)]xdaylight(35",Mar2\)  F^a(h,atb,vis) 

cos5[0(to,0<7y)]  x  daylightjlal, day)  ,  El  (h, alb , vis) 

cos5  [0(35°,  A  for  2 1)]  x  daylight  ( 35” ,  Xfo^ti)  ^ ( E?  (b,  alb ,  vis) } 


cos4[0(/o/,  day))  x  daylighiilat ,  day) 
cos4[0(35“,  Afor  21)]  x  daylighi(i5\Mar2\) 


for  v  =  E'f  ( h,  alb ,  \ is) 


(13) 


where  e  is  the  solar  zenith  angle  at  noon  (shown  in  eq  3  for  hour  =  1 2)  as  a  function  of  latitude, 
and  daylight  is  the  length  of  daylight  in  hours  (eq  3).  For  example,  the  total  daily  action-energy 
for  45°latitude  on  1  August  for  1976  U.S.  standard  atmosphere  with  rural  aerosol  model  at 
spring-summer  (Table  1)  is  given  by  eq  13  as 

El (h, alb ,  vis ,45" ,  Junl 1)  =  1 .246  x  (55.7  +  17.6  h-  0.923  h2  +  7 1 .2  alb  +  28.7  alb 2  +  0. 1 48  vis  -  0.000888  vis2 ) 

where  the  magnitude  scaling  constant  of  (1.246)  cos  l^45  » Alls  2 hi x  daylight (45  ,  Aug  21)  =  {  246 

cos5[0(35".  Mar  21)]  x  daylighl(35°,Mar2\) 

is  computed  from  eq  3  for  21  August  (233rd  day  of  the  year)  at  45°  latitude  where  0  =  32.63°  (at 
noon)  and  daylight  =  13.69  hours ;  and  for  the  reference  2 1  March  (80th  day  of  the  year)  at  35° 

latitude  where  0  =  35. 1 1°  (noon). 


We  emphasize  that  this  is  only  a  magnitude  scaling,  and  it  does  not  reshape  the 

dependence  (vertical  profile)  of  the  variable  (i.e.,  the  ratio,  - y(h,aib,  vis,  latitude,  day ) - ^  js  a 

y(h,alb,  vis ,1  at  -  35 ° ,  day  -  Mar  21) 

constant  that  is  not  a  function  of  h,  alb,  and  vis).  This  is  a  deficiency  in  the  location/time  scaling, 
but  it  produces  reasonable  results  as  is  evident  in  Table  5.  In  Table  5,  we  show  the  error  in 
latitude/time  scaling  for  1976  U.S.  standard  atmosphere  and  mid-latitude  summer  atmosphere 
with  rural  and  urban  aerosols  (noted  as  atmospheres  A,  B,  C,  and  D)  when  we  scale  the  reference 
atmosphere  (21  March  and  35°  latitude)  to  latitudes  25,  35,  and  45°  for  four  dates  (1  February, 

1  May,  1  August,  and  1  November).  The  scaling  results  for  these  six  quantities  are  very 
satisfactory  for  some  days/latitudes  (<5%  error),  but  for  the  other  days/latitudes,  the  errors  are  as 
high  as  50%  (e.g.,  -50%  for  Efrand  -40%  for  El  at  February  and  November  for  45°  latitude). 
The  average  of  all  288  entries  in  Table  5  is  -10%,  whereas  the  error  is  <5%  for  130  entries. 
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Thus,  the  simple  approach  of  scaling  to  predict  the  solar  UV  irradiance  received  by  airborne  viral 
pathogens  in  the  atmosphere  with  <50%  error  should  be  useful  for  theoretical  studies  of  UV 
inactivation  of  viruses  in  the  atmosphere,  where  virus  survival  is  estimated  in  logarithmic  scale 
with  large  uncertainty  ( ±0.51ogl0). 


5.  CONCLUSIONS 

c 

Inactivation  of  viruses  suspended  in  clear-sky  atmosphere  due  to  direct  and 
diffuse  solar  UV  was  previously  (/)  calculated  with  MODTRAN  plane  irradianccs.  However,  / 

MODTRAN  is  a  complex  coded  software  program,  and  computing  daily  energies  is  time 
consuming.  Therefore,  we  developed  a  simple  nonlinear  regression  model  that  can  produce 
various  quantities  (eqs  6  and  7)  within  a  reasonable  accuracy  (  ~-5%  of  the  values  produced  by 
MODTRAN)  that  arc  useful  to  theoretical  and  experimental  UV  inactivation  studies. 

Development  of  a  regression  model  also  provided  an  insight  into  the  dependence  of  key 
parameters  (altitude,  albedo,  and  visibility),  whereas  MODTRAN  (a  numerical  code)  “only” 
provided  numerical  (though  extremely  accurate)  output.  To  the  best  of  our  knowledge,  this  was 
the  first  time  that  a  comprehensive  regression  model  was  demonstrated  on  MOD  I  RAN 
atmospheric  radiances.  With  the  spectral  shape  of  the  action  spectrum  and  the  narrow 
wavelength  range  of  280-320  nm,  we  were  pleasantly  surprised  that  a  reasonable  accuracy  in  the 
regression  model  was  accomplished  for  our  application  (Tables  1—4,  last  two  columns). 

We  chose  nine  quantities  (listed  in  cq  6  and  shown  in  Tables  1—4)  for  the 
regression  model  (normalized  by  ambient  radiation  fields).  These  quantities  may  be  useful  in 
experimental  work  when  one  wants  to  rescale  measured  ambient  data,  and  the  absolute  quantities 
may  be  useful  for  theoretical  studies.  The  nine  quantities  were  ETaXEf  XEJ"  (total,  diffuse  and 
direct  daily  action-energies),  ET0IFlxx,  E?  /F^,  ETa  /  F^x  (daily  action-energies  that  are 
normalized  by  atmospheric  ambient  tlux  measured  at  noon),  Edi/  /  ,  EJJr  /  ETa  ,  and 

F^x  /  F2XX  (fraction  diffuse  and  direct  daily  action-energy,  and  fraction  of  diffuse  ambient  flux  at 
noon).  From  these  nine  quantities,  we  derived  four  additional  ambient  quantities  (eq  7): 

(total,  diffuse,  and  direct  ambient  fluxes  at  noon),  and  EJir (ambient  direct  daily 

energy). 


We  developed  a  second  order  regression  model  (to  capture  the  nonlinearity  of  the 
radiation  quantities  (the  Figure)  in  the  form  of 

\ih, alb, vis)  =  h,  +  u-,/j  +  w\h2  +  w4alb  +  w,alb2  +  wtvis  +  w^vis2  (eq  9,  Tables  1—4)  for  the  nine  quantities 
for  atmospheric  models,  1976  U.S.  standard  and  mid-latitude  atmospheres  for  two  seasons 
(spring- summer  and  fall-winter)  and  four  dates  (equinox  and  solstice).  The  regression  model 
reproduced  MODTRAN  radiation  quantities  with  acceptable  accuracy  (mean  accuracy  of  ~2%); 
the  worst  accuracy  was  ~6%  for  radiation  fields  E?r  and  Fddx  / eI  on  21  December.  For  all  other 
quantities,  the  error  was  <5%.  This  accuracy  is  excellent,  compared  to  microbiological 
laboratory  experiments  where  the  survival  level  of  viruses  to  UV  radiation  was  measured  in  log- 
scale  units,  usually  within  an  accuracy  no  better  than  ±0.5  logio  units.  With  our  simple  regression 
model,  we  can  compute  (cq  10)  the  maximum  number  of  radiation-days,  needed  for 
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inactivating  virus  X  at  height  h  in  a  given  atmospheric  model,  visibility  condition,  and  ground 
albedo. 


The  accuracy  of  our  regression  model  was  satisfactory  (Tables  1 — 4),  but  its  utility 
was  limited  to  a  specific  latitude  (35°)  and  day  of  the  year  (equinox  and  solstice  dates).  One 
solution  was  to  recompute  regression  coefficients  (eq  12)  for  the  desired  latitude/date,  a  process 
that  was  time  consuming  and  involved  a  vast  number  of  MODTRAN  runs.  As  an  alternative,  we 
attempted  to  scale  (see  eq  13)  the  regression  model  from  a  reference  latitude  and  date  (chosen  as 
35°  latitude  for  21  March)  to  an  arbitrary  latitude  and  day  of  the  year.  We  found  a  simple  scaling 
as  a  function  of  zenith  angle  and  length  of  daylight  for  the  six  quantities,  El  I  F^x ,  Ef  /  , 

F^%  /F,lax,  Era ,  E*{ ,  E*r .  The  accuracy  of  our  simple  location/time-scaling  approach  was  less 
viable  (see  Table  5)  than  the  excellent  accuracy  of  the  regression  model  (see  Tables  1—4),  but  it 
was  still  within  the  general  accuracy  attained  in  the  virus  inactivation  studies. 

The  approach,  equations,  and  procedure  presented  in  this  study  should  assist  in 
estimating  the  time  during  which  viral  and  other  pathogens  remain  infectious  after  accidental  or 
intentional  release  in  the  atmosphere.  The  demonstrated  success  of  the  regression  model  should 
encourage  other  researchers  to  seek  simple  regression  models  for  complex  MODTRAN 
computations  in  a  variety  of  additional  applications. 


Table  !.  1976  U.S.  Standard  Atmosphere  with  Rural  Aerosols  at  Latitude  35°  at  Four  Equinox 
and  Solstiee  Dates 


7  1  2 

Regression  Model:  y(h,alb,vis)  =  »v,  +  w^h  +  +  wAaib  +  w^alb*  +  w6v/.y  +  w^vis 

0  <>  h  <>  5 (Aw)  0  £  albedo  <>  0. 1 5  1 5  Aw  <  visibility  £  1 00 Aw 


y 

Day 

w 7 

\V2 

M’J 

HV 

H’J 

U6 

H’7 

Mean 
(error  | 
(%) 

Max 

1  error) 

(%) 

with 

P  >  0.9 

Et 

^  a 

Mar  21 

2.6IE+01 

9.79E-01 

-5.20E-02 

-1.87E+00 

6.60E-02 

1  17E-02 

-7.08E-05 

0.18 

0.37 

Jun  21 

4.0IE+0I 

1 .66E+00 

-9.59E-02 

-3.36E+00 

1.69E-01 

1  90E-02 

-I.I5E-04 

0.22 

0.47 

ft 

1  max 

Sep  21 

2.68E+0I 

1 03E+00 

-5.73E-02 

-1.95E^00 

7.88E-02 

1.18E-02 

-7.I4E-05 

0.19 

0.39 

Dec  21 

1.14E+01 

2.20E-01 

-4.26E-03 

-4. 85  E-01 

-1.37E-02 

2.49E-03 

-1.47E-05 

0.09 

0.20 

co 

Mar  21 

2.77E+01 

1.10E+00 

-8.48E-02 

-7.60E+00 

6.66E+00 

3.78E-02 

-2.21  E-04 

0.41 

0.94 

Jun  21 

4.30E+01 

1.80E+00 

-1.56E-0I 

-1.44E+0I 

1.46E+01 

5.32E-02 

-3.15E-04 

0.43 

0.94 

ctJif 

rmax 

Sep  21 

2.86E+0I 

1.15E+00 

-9.60E-02 

-8.34E+00 

7.69E+00 

3.88E-02 

-2.28E-04 

0.43 

0.95 

Dec  21 

1.18E+0I 

3.43  E-01 

-1 .64  E-02 

-1.91E+00 

1.22E+00 

1 .59E-02 

-9.00E-05 

0.33 

0.78 

c<hr 

"max 

Mar  21 

8.58E-03 

8.38E-04 

-2.58E-05 

-1.43E-02 

9.13E-03 

I.52E-04 

-8.8  IE-07 

4.57 

10.05 

Jun  21 

7.32E-03 

3.73E-04 

-7.43E-06 

-I.10E-02 

7.I2E-03 

9.93E-05 

-5.82E-07 

3.96 

9.14 

eI 

Sep  21 

9.08E-03 

8.80E-04 

-4.50E-05 

-1.47E-02 

9.37E-03 

1  49E-04 

-8.70E-07 

4.50 

9.81 

Dec  2 1 

7.65E-03 

3.42E-03 

- 1 .46E-04 

-1 .93E-02 

1.18E-02 

2.83E-04 

-1  60E-06 

6.38 

13.85 

El 

Mar  21 

I.2IE+00 

6.47E-02 

3.87E-03 

-7.51  E-01 

7.28E-0I 

6.96E-03 

-3.8IE-05 

1.55 

3.65 

Jun  21 

1.32E+00 

7.17E-02 

4.66E-03 

-I.16E+00 

1 .28E+00 

8.92E-03 

-4.89E-05 

1.87 

4.30 

L"  dl  f 

c-a 

Sep  21 

1.23E+00 

7.61  E-02 

2.48E-03 

-8.47E-0I 

8.46E-01 

7.41  E-03 

-4.09E-05 

1.68 

3.86 

Dec  21 

1 .05E+00 

4.90E-02 

2.1  IE-03 

-2.46E-0I 

1.79E-01 

3.30E-03 

-1.79E-05 

0.86 

2.05 

EJir 

Mar  21 

1.34E+02 

2.24E+01 

-5.75E-0I 

-2.47E+02 

1.56E+02 

2.91  E+00 

-1.67E-02 

1.55 

3.65 

Jun  21 

1.38E+02 

I.48E+01 

-3.80E-0I 

-2.28E+02 

1  46E+02 

2.33E+00 

-1.35E-02 

1.87 

4.30 

El 

Sep  21 

I.44E+02 

2.41E+01 

-I.05E+00 

-2.58E+02 

1.63E+02 

2.90E-*-00 

-1.68E-02 

1.68 

3.86 

Dec  21 

8.42  E+01 

5.35E+01 

-I.45E+00 

-2.50E+02 

1.50E+02 

3.91  E+00 

-2.18E-02 

0.86 

2.05 

pdif 

max 

Mar  21 

7.77E-01 

-3.48E-02 

9.I9E-04 

4.14E-01 

-2.65E-01 

-4.I7E-03 

2.42E-05 

2.45 

5.75 

Jun  21 

7.08E-01 

-3. 15E-02 

9.04  E-04 

4.95  E-01 

-3.22E-01 

-4.24E-03 

2.48E-05 

2.81 

6.56 

ft 

1  max 

Sep  21 

7.58E-0 1 

-3.77E-02 

1.60E-03 

4.36E-0I 

-2.80E-0I 

-4.20E-03 

2.46E-05 

2.63 

6.17 

Dec  21 

9.14E-01 

-4.32  E-02 

1.31E-03 

2.35E-0I 

-1.44E-01 

-3.31  E-03 

1.87E-05 

1.58 

3.74 

El 

Mar  21 

5.57E+01 

1.76E+01 

-9.23E-01 

7.I2E+0I 

2.87E+01 

I.48E-0I 

-8.88E-04 

1.74 

3.58 

Jun  21 

I.22E+02 

3.47E+01 

-I.92E+00 

E53E+02 

6.02E+01 

2.90E-01 

-1.74E-03 

1.57 

3.24 

Sep  21 

5.87E+01 

1.84E+0I 

-1.01E+00 

7.45E+0I 

2.98E+01 

1.46E-01 

-8.79E-04 

1.72 

3.54 

Dec  21 

9.1 1E+00 

3.29E+O0 

-1.49E-01 

1.20E+01 

4.99E+00 

2.67E-02 

-1 .60E-04 

1.96 

4.00 

Ef 

Mar  2 1 

4.73E+01 

1  00E+01 

-1.02E+00 

7.12E+0I 

2.87E+0I 

-1.53  E-01 

8.16E-04 

1.23 

2.52 

Jun  21 

9.49E+01 

1.72E+01 

-1.87E+00 

1.53E+02 

6.02E+01 

-4.05E-01 

2.24E-03 

1.22 

2.49 

Sep  21 

4.88E+01 

9.78E+00 

-1.00E+00 

7.45E+01 

2.98E+0I 

-1.73  E-01 

9. 42  E-04 

1.22 

2.50 

Dec  21 

8.8IE+00 

2.40E+00 

-2.07E-0I 

I.20E+01 

4.99E+00 

-5.49E-03 

1.65E-05 

1.56 

3.20 

r'tJlf 

ba 

Mar  21 

8.36E+00 

7.60E+00 

I.02E-01 

-2.51  E- 1 3 

I.46E-I2 

3.01  E-01 

-1.70E-03 

4.82 

10.26 

Jun  21 

2.67E+01 

1.75E+OI 

-5.39E-02 

-6.76E-13 

3.75E-I2 

6.95E-01 

-3.98E-03 

4.32 

9.37 

Sep  21 

9.90E+00 

8.65E+00 

-1.09E-02 

-1.85E-13 

1.05E-I2 

3.18E-01 

-1.82  E-03 

4.89 

10.31 

Dec  21 

3.04  E-01 

8.95E-0I 

5.74E-02 

-3.53E-14 

I.92E-I3 

3.22E-02 

-I.76E-04 

6.37 

15.25 

Altitude  h  and  visibility  vis  are  in  kilometers.  Albedo  alb  is  in  fractions  (0-1  ).This  is  a  rural  aerosol  model 
(MODTRAN  parameter  ihaze).  21  March  and  2 1  June  are  with  the  spring-summer  seasonal  model  (MODTRAN 


parameter  iseasn)\  21  September  and  21  December  are  with  fall-winter  seasonal  model.  The  mean  absolute  error  in 
the  regression  model  and  the  maximum  absolute  error  (for  probability  >0.9)  were  computed  for  1760  values  /j,  alb, 
and  vis.  The  Table  shows  excellenl  fil  of  the  regression  model  lo  the  complex  MODTRAN  radiation  fields  and  2  1% 
mean  error  (for  all  36  entries).  The  worst  accuracy  (-6%)  is  for  radiation  fields  Efr  and  F^x  /  E*  at  2 1  December. 

For  all  other  radiation  fields,  lhe  mean  absolute  error  is  less  lhan  5%  (for  mosl  enlries  the  error  is  <2%,  mean  of 
1.85%  for  all  other  34  enlries). 
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Table  2.  1976  U  S.  Standard  Atmosphere  with  Urban  Aerosols  at  Latitude  35°  at  Four 
Equinox  and  Solstice  Dates _ 

Regression  Model:  y(h,alb,vis)  =  Hj  +  Wjh  +  w^h~  +  w^alb  +  w$alb2  +  w^vis  +  Wjvis‘ 

0  <  h  <  5 (km)  0  <,  albedo  <  0. 1 5  1 5  km  <  visibility  <  1 00  km 


y 

Day 

Wj 

W2 

W3 

W}4 

Vt;5 

Vi’6 

W; 

Mean 

\error\ 

(%) 

Max 
|  error | 
(%) 
with 

P>  0.9 

Et 

^ a 

FJ 

1  max 

Mar  21 

2.60E+01 

1.03E+00 

-5.51  E-02 

-1.62E+00 

1.49E-01 

1  27E-02 

-7.72E-05 

0.25 

0.57 

Jun  21 

3.98E+01 

1.74E+00 

-1.01E-01 

-2.90E+00 

3.04E-01 

2.27E-02 

-1.39E-04 

0.30 

0.70 

Sep  2! 

2.67E+01 

1.08E+00 

-6.05E-02 

-1.68E+00 

1.63E-01 

1 .32E-02 

-8.06E-05 

0.26 

0.60 

Dec  21 

1.15E+01 

2.25E-01 

-4. 1 7E-03 

-4.13E-01 

1.46E-02 

1.21E-04 

-3.34E-08 

0.08 

0.16 

Ed>f 

c  a 

rJ'f 

max 

Mar  21 

2.83E+01 

1.11E+00 

-8.21E-02 

-7.57E+00 

7.07E+00 

2.73E-02 

-1.62E-04 

0.48 

1.10 

Jun  21 

4.37  E+01 

1.84E+00 

-1.53E-01 

-1.43E+01 

1.53E+01 

3.92  E-02 

-2.37E-04 

0.50 

1.16 

Sep  21 

2.92E+01 

1.17E+00 

-9.34E-02 

-8.35E+00 

8.18E+00 

2.76E-02 

-1.66E-04 

0.50 

1.14 

Dec  21 

1.22E+01 

3.37E-01 

-1.5  IE-02 

-1.90E+00 

1.32E+00 

9.68E-03 

-5.50E-05 

0.34 

0.78 

r 'dit 
~max 

et 

Mar  21 

1 . 1 6  E-02 

4.90E-04 

7.72E-06 

-1.44E-02 

9.30E-03 

9.81E-05 

-5.71E-07 

3.27 

7.21 

Jun  21 

9.47E-03 

1.07E-04 

1.95E-05 

-1.10E-02 

7.20E-03 

5.93E-05 

-3.46E-07 

2.73 

6.31 

Sep  21 

1.2  2  E-02 

5.21E-04 

-1.03E-05 

-1.48E-02 

9.55E-03 

9.26E-05 

-5.45E-07 

3.19 

7.15 

Dec  21 

1.16E-02 

3.1  IE-03 

-1.28E-04 

-1.96E-02 

1.22E-02 

2.20E-04 

-1.27E-06 

5.20 

11.18 

ETa 

Ef 

Mar  21 

1.32E+00 

5.63E-02 

4.71E-03 

-7.97E-01 

7.91E-01 

5.23E-03 

-2.91E-05 

1.48 

3.37 

Jun  21 

1.48E+00 

5.80E-02 

6. 1 0E-03 

-1.25E+00 

1.40E+00 

6.47E-03 

-3.60E-05 

1.74 

3.90 

Sep  21 

1.36E+00 

6. 73  E-02 

3.30E-03 

-9. 04  E-01 

9.25E-01 

5.42E-03 

-3.05  E-05 

1.60 

3.57 

Dec  21 

1.09E+00 

4.84E-02 

2. 12E-03 

-2.54E-01 

1.89E-01 

2.72E-03 

-1.51E-05 

0.92 

2.16 

EJ" 

El 

Mar  21 

1.83F-02 

1.70E+01 

-8.29E-02 

-2.49E+02 

1.59E+02 

2.05E+00 

-1.18E-02 

1.48 

3.37 

Jun  21 

1.81E+02 

9.79E+00 

1.08E-01 

-2.29E+02 

1 .48E+02 

1.55E+00 

-8.99E-03 

1.74 

3.90 

Sep  21 

1 .96E+02 

1 . 84E+0 1 

-5.39E-01 

-2.60E+02 

1.66E+02 

1.98E+00 

-1.16E-02 

1.60 

3.57 

Dec  21 

1.32E+02 

5. 04  E+01 

-1.32E+00 

-2.53E+02 

1.55E+02 

3.16E+00 

-1.80E-02 

0.92 

2.16 

pdif 

1  max 

ft 

1  max 

Mar  21 

7.00E-01 

-2.71E-02 

1.76E-04 

4.14E-01 

-2.67E-01 

-2.76E-03 

1.60E-05 

2.12 

4.86 

Jun  21 

6.22E-01 

-2.26E-02 

-7.66E-06 

4.91E-01 

-3.22E-01 

-2.63E-03 

1.54E-05 

2.35 

5.38 

Sep  21 

6.75  E-01 

-2.96E-02 

8.23E-04 

4.36E-01 

-2. 82  E-01 

-2.68E-03 

1.58E-05 

2.27 

5.24 

Dec  21 

8.67E-01 

^t.02E-02 

1. 12E-03 

2.38E-01 

-1.48E-01 

-2.55E-03 

1.47E-05 

1.54 

3.61 

et 

c a 

Mar  21 

4.09E+0 1 

1.81E+01 

-9.02  E-01 

6.21  E+01 

2.09E+01 

4.66  E-01 

-2.80E-03 

2.46 

4  81 

Jun  21 

9.25E+01 

3.59E+01 

-1.93E+00 

1.34E+02 

4.42E+01 

9.18E-01 

-5.54E-03 

2.26 

4.37 

Sep  21 

4.33E+01 

1.90E+0I 

-9.96E-0 1 

6.50E+01 

2.17E+01 

4.78E-01 

-2.88E-03 

2.43 

4.73 

Dec  2 1 

6.33E+00 

3.34E+00 

-1.40E-01 

1.03E+01 

3.58E+00 

8. 60  E-02 

-5.13E-04 

2.69 

5.29 

edif 

Mar  21 

3.20E+01 

1.07E+01 

-1.02E+00 

6.21  E+01 

2. 09 E+01 

1.72E-01 

-1.13E-03 

1.89 

3.73 

Jun  21 

6.46E+01 

1.88E+01 

-1.92E+00 

I.34E+02 

4.42E+01 

2.42E-01 

-1.65E-03 

1.72 

3.52 

Sep  21 

3.28E+01 

1.05  E+01 

-1.01E+00 

6.50E+01 

2.17E+01 

1.68E-01 

-1.10E-03 

1.82 

3.55 

Dec  21 

5.97E+00 

2.46E+00 

-1.99E-01 

1.03E+01 

3.58E+00 

5. 43  E-02 

-3.39E-04 

2.33 

4.55 

r-dir 

ba 

Mar  21 

8  87E+00 

7.45E+00 

1.19E-01 

-2.40E-13 

1.42E-12 

2.94E-01 

-1.67E-03 

4.58 

9.90 

Jun  21 

2.79E+01 

1.71E+01 

-1.01E-02 

-4.93E-13 

2.95E-12 

6.76E-01 

-3.88E-03 

4.11 

8.91 

Sep  21 

1.05E+01 

8.49E+00 

8.95E-03 

-2.69E-13 

1.36E-12 

3.10E-01 

-1.78E-03 

4.64 

9.96 

Dec  21 

3.52E-01 

8.81E-01 

5.89E-02 

-1.46E-14 

1.05E-13 

3.17E-02 

-1.75E-04 

6.03 

14.56 

Same  as  Table  1 ,  but  for  1976  U.S.  standard  atmosphere  with  urban  aerosols.  The  Table  shows  excellent  fit  of  the 
regression  model  and  2%  mean  error  (all  36  entries).  The  worst  accuracy  (~6%)  is  for  radiation  fields  Ef  and 

/■max  / Etq  at  21  December.  For  all  other  radiation  fields,  the  mean  absolute  error  <5%  (mean  of  1 .84%  for  all  other  34 
entries). 


27 


Table  3.  Mid-Latitude  Atmosphere  with  Rural  Aerosols  at  Latitude  35°  at  Four  Equinox  and 
Solstice  Dates 


Regression  Model:  y(h9alb>vis)  =  m(  +  wyh  +  w^h2  +  \\Aa!h  +  w$albm  +  vv6W.v  +  wjvis2 
0  <h<  5 {bn)  0  <  albedo  <0.15  1 5km  <1  visibility  <i  1 00km 


y 

Day 

W/ 

\V3 

ny 

U’5 

W6 

U’7 

Mean 

error 

(%) 

Max 

error 

<%) 

with 

P>  0.9 

Et 

r'a 

ft 

rtTMUC 

Mar  21 

2.69E+01 

I.01E+00 

-4.04E-02 

-2.21E400 

3.51  E-02 

1.27E-02 

-7.62E-05 

0.19 

0.38 

Jun.  21 

4.14E+01 

I.70E+00 

-7.65E-02 

-3.94E400 

1.31  E-01 

2.03E-02 

-1.23E-04 

023 

0.47 

Sep  21 

2.33E+01 

8.62E-01 

-5.25E-02 

-1.61E400 

6.65E-02 

9.47E-03 

-5.74E-05 

0.19 

0.38 

Dec  21 

1.0IE+01 

1.79E-01 

-4.13E-03 

-4.02E-01 

-1.33E-02 

1.99E-03 

-1.18E-05 

0.09 

0.19 

pdif 

‘'a 

‘  max 

Mar  21 

2.86E+01 

1.1  IE+00 

-7.91  E-02 

-8.09E400 

6.85  E400 

3.79E-02 

-2.22E-04 

0.40 

0.88 

Jun  21 

4.43E+01 

1.79E400 

-1.47E-01 

-I.52E+01 

1.490+01 

5.17E-02 

-3.07E-04 

0.41 

0.88 

Sep  21 

2.48E+0I 

9.69E-01 

-8.98E-02 

-6.98E400 

6.47E400 

3.20E-02 

-I.88E-04 

0.42 

0.94 

Dec  21 

1.05E+01 

2.93E-01 

-1.51  E-02 

-1.66E+00 

1.07E400 

I.40E-02 

-7.88E-05 

0.34 

0.78 

rp  <Jtr 

''max 

El 

Mar.  21 

8.32E-03 

8.27E-04 

-3.56E-05 

-1.37E-02 

8.77E-03 

1 .48E-04 

-8.57E-07 

4.61 

10.21 

Jun.  21 

7.10E-03 

3.75E-04 

-1.51E-05 

-1.05E-02 

6.83  E-03 

9.67E-05 

-5.67E-07 

4.00 

9.28 

Sep  21 

1.05E-02 

1.09E-03 

-5.07E-05 

-I.70E-02 

1.08E-02 

1.73E-04 

-I.02E-06 

4.53 

9.80 

Dec  21 

8.63E-03 

4.04E-03 

-1.65E-04 

-2.19E-02 

1.33E-02 

3.23E-04 

- 1 .83  E-06 

6.38 

13.83 

El 

F-f 

Mar  21 

I.21E+00 

6.69E-02 

3.96E-03 

-7.61  E-01 

7.40E-01 

7.13E-03 

-3.90E-05 

1.58 

3.71 

Jun  21 

1.32E+00 

7.54E-02 

4.80E-03 

-1.19E+00 

1.31  E+00 

9.20E-03 

-5.04E-05 

1.91 

4.39 

Sep  21 

1.24E+00 

7.90E-02 

3.23E-03 

-8.70E-01 

8.70E-0I 

7.59E-03 

-4.19E-05 

1.70 

3.90 

Dec  21 

I.05E+00 

5.04E-02 

2.35E-03 

-2.48E-01 

1.80E-0I 

3.33E-03 

-1.81E-05 

0.86 

2.06 

EJ" 

Mar  21 

1.30E+02 

2.I9E+01 

-7.66E-0I 

-2.37E+02 

1.50E+02 

2.83E400 

-1  63  E-02 

158 

3.71 

Jun  21 

1.34E+02 

1.46E+01 

-5.51E-01 

-2.18E+02 

1.40E402 

2.27E4  00 

-1.32  E-02 

1.91 

4.39 

El 

Sep  21 

1.66E+02 

2  91E+01 

-I.I8E+00 

-2.97E+02 

1.87E+02 

3.37E+00 

-1.95  E-02 

1.70 

3.90 

Dec  21 

9.47E+01 

6.26E+01 

-  I.57E+00 

-2.83E+02 

1.69E+02 

4.43E+00 

-2.48E-02 

0.86 

2.06 

1  max 

ft 

1  max 

Mar  21 

7.77E-01 

-3.52E-02 

1.03E-03 

4.14E-01 

-2.66E-0I 

-4.20E-03 

2.43E-05 

2.47 

5.80 

Jun  21 

7.07E-01 

-3.2  IE-02 

1.03E-03 

4.95E-01 

-3.23E0I 

-4.26E-03 

2.49E-05 

2.84 

6.64 

Sep  21 

7.57E-01 

-3.91  E-02 

1.61E-03 

4.37E-01 

-2.7QE-01 

-4.24E-03 

2.48E-05 

2.70 

6.31 

Dec  21 

9.14E-0I 

-4.47E-02 

1.32E-03 

2.36E-01 

-1.43  E-01 

-3.34E-03 

1.89E-05 

1.61 

3.79 

et 

L'a 

Mar  21 

5.82E+01 

I.79E+01 

-8.45E-0I 

7.33E+01 

2.83E401 

I.52E-01 

-9.I2E-04 

1.72 

3.54 

Jun  21 

1.27E+02 

3.52E+01 

-I.76E+00 

1.57E402 

5.90E+01 

2.97E-01 

-I.78E-03 

1.55 

3.21 

Sep  21 

4.79E+0I 

1 .51E+01 

-8.86E-01 

6.03E401 

2.41E40I 

I.11E-01 

-6.65E-04 

1 .65 

3.39 

Dec  21 

7.44  E+00 

2.71E+00 

-1.3  IE-01 

9.69E400 

4.03E+00 

2.05E-02 

-1.22E-04 

1.87 

3.87 

p-dif 

Mar  21 

4.93E-*01 

I.01E+01 

-9.94E-01 

7.33E+01 

2.83E401 

-1.65  E-01 

8.86E-04 

1.22 

2.51 

Jun  21 

9.87E+0I 

I.70E+0I 

-1.80E+00 

1.57E+02 

5.90E401 

-4.360-01 

2.42E-03 

1.25 

2.56 

Sep  21 

3.97E+0I 

7.89E+00 

-8.65E-01 

6.03E^01 

2.4  IE  401 

-1.50E-01 

8.26E-04 

1.24 

2.54 

Dec  21 

7.19E+00 

1 .96E4  00 

-1.78E-01 

9.69E400 

4.03E400 

-5.79E-03 

2.19E-05 

1.50 

3.08 

edir 

Mar  21 

8.85E+00 

7.89E400 

1.49E-01 

-2.33E-13 

1.33E-12 

3.18E-01 

-1.80E-03 

4.81 

10.36 

Jun  21 

2.82E-I0I 

I.82E+OI 

4.37E-02 

-7.13E-13 

3.82E-12 

7.34E-01 

-4.21  E-03 

4.31 

9.40 

Sep  21 

8.13E+00 

7.23E400 

-2.08E-02 

-2  76E-13 

1.38E-I2 

2.60E-01 

-1.49E-03 

4.84 

10.27 

Dec  21 

2.47E-01 

7.48E-01 

4.64E-02 

-1.08E-14 

5.29E-14 

2.62E-02 

-1  44E-04 

6.28 

14.99 

Same  as  Table  1  but  for  mid-latitude  atmosphere  with  rural  aerosols.  2 1  March  and  21  June  are  for  mid-latitude- 
summer  atmosphere  at  spring-summer  season  model.  2 1  September  and  2 1  December  are  for  mid-latitude-winler 
atmosphere  at  spring-summer  season.  The  Table  shows  excellent  fit  of  lhe  regression  model  and  2.1%  mean  error 


(considering  all  36  entries).  The  worst  accuracy  (-6%)  is  for  radiation  fields  E%r  and  F^x  IF*  at  2 1  December.  For 
all  olher  radialion  fields,  lhe  mean  absolute  error  is  <5%  (a  mean  of  1 .85%  for  all  other  34  entries). 
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Table  4.  Mid-Latitude  Atmosphere  with  Urban  Aerosols  at  Latitude  35°  at  Four  Equinox  and 
Solstice  Dates 


Regression  Model:  y(h,alb,vis)  -  *'\  +  +  u^/r  +  wAalb  +  w$a!b2  +  w^vis  +  h^vw2 

0  <,h£  5(km)  0  ^  albedo  ^  0. 1 5  1 5  km  £  visibility  £  1 00 Am 


y 

Day 

w, 

H’2 

H’i 

MV 

M’J 

M  ’6 

H7 

Mean 

' error 

(%) 

Max 
l  error] 
(%) 
with 

P>  0.9 

ET. 

Mar  21 

2.69E+01 

1.06E+00 

-4.30E-02 

-1.93E+00 

1.32E-01 

I.20E-02 

-7.30E-05 

0.24 

0.55 

Jun  21 

4.13E+1 

1.78E+0 

-8.06E-2 

-3.43E+0 

2.88E-1 

2.13E-2 

-1.31E-4 

0.29 

0.67 

eL 

Sep  21 

2.33E+01 

8.62  E-01 

-5.25  E-02 

-1.61E+00 

6.65E-02 

9.47E-03 

-5.74E-05 

0.19 

0.38 

Dec  21 

1.02E+0I 

1.82E-01 

-4.02E-03 

-3.49  E-01 

8.83E-03 

1.53E-04 

-3.34E-07 

0.08 

0.16 

E« 

Mar  21 

2.92E+01 

1.12E+00 

-7.67E-02 

-8.01  E+00 

7.24E+00 

2.63E-02 

-1.57E-04 

0.46 

1.07 

Jun  21 

4.5 1E+01 

1.83E+00 

-1.43E-01 

-I.50E+0I 

1.56E+01 

3.58E-02 

-2.19E-04 

0.48 

1.1 1 

pdf 

1  max 

Sep  21 

2.48E+0t 

9.69E-01 

-8.98E-02 

-6.98E+00 

6.47E+00 

3.20E-02 

-t.88E-04 

0.42 

0.94 

Dec  21 

1.07E+01 

2.90E-01 

-1.43E-02 

-1.66E+O0 

1.15E+00 

9.15E-03 

-5.19E-05 

0.35 

0.80 

cdir 
*  max 

Mar  21 

1.09E-02 

5.27E-04 

-7.79E-06 

-1.38E-02 

8.93E-03 

1.01E-04 

-5.86E-07 

3.43 

7.59 

Jun  21 

8.99E-03 

1 .45E-04 

7.54E-06 

-1.06E-02 

6.91E-03 

6.16E-05 

-3.60E-07 

2.87 

6.59 

El 

Sep  21 

I.05E-02 

1.09E-03 

-5.07E-05 

-I.70E-02 

I.08E-02 

I.73E-04 

-I.02E-06 

4.53 

9.80 

Dec  21 

I.25E-02 

3.78E-03 

-1.53E-04 

-2.22E-02 

I.37E-02 

2.60E-04 

-1.50E-06 

5.40 

11.63 

El 

Mar  21 

1.32E-MX) 

5.98E-02 

4.62E-03 

-8.04E-01 

7.99E-01 

5.55E-03 

-3.08E-05 

1.53 

3.48 

Jun  21 

1.47E+00 

6.35E-02 

5.97E-03 

-1.26E+00 

1 .43E+00 

6.96E-03 

-3.87E-05 

1.80 

4.04 

Ef 

Sep  21 

1.24E+00 

7.90E-02 

3.23E-03 

-8.70E-01 

8.70E-01 

7.59E-03 

-4.19E-05 

1.70 

3.90 

Dec  21 

1.09E+00 

5.04E-02 

2.30E-03 

-2.55E-01 

1.89E-01 

2.81  E-03 

-1.56E-05 

0.94 

2.20 

Edir 

Mar  2 1 

1.73E+02 

1.73E+01 

-3.62E-0I 

-2.39E+02 

1.53E+02 

2.08E+00 

-1.20E-02 

1.53 

3.48 

Jun  21 

1.72E+02 

1 .03E+01 

-1.43E-01 

-2.19E+02 

1.42E+02 

1.59E+00 

-9.20E-03 

1.80 

4.04 

El 

Sep  21 

1.66E+02 

2.91E+01 

-1.1 8E+00 

-2.97E+02 

1.87E+02 

3.37E+00 

-1.95E-02 

1.70 

3.90 

Dec  21 

1  42E+02 

6.02E+01 

-1.54E+00 

-2.86E+02 

1.74E+02 

3.70E+00 

-2.10E-02 

0.94 

2.20 

1  max 

Mar  21 

7.06E-01 

-2.84F.-02 

4.03E-04 

4.14E-01 

-2.68E-0I 

-2.90E-03 

1 .69E-05 

2.19 

5.04 

Jun  21 

6.28E-01 

-2.41E-02 

2.44E-04 

4.92E-01 

-3.23E-01 

-2.79E-03 

I.63E-05 

2.44 

5.56 

ft 

1  max 

Sep  21 

7.57E-01 

-3.91E-02 

1.61E-03 

4.37E-01 

-2.79E-01 

-4.24E-03 

2.48E-05 

2.70 

6.31 

Dec  21 

8.72E-01 

-4.24E-02 

1.23E-03 

2.38E-01 

-1.47E-01 

-2.67E-03 

1.53E-05 

1.60 

3.74 

et 

Mar  21 

4.40E+01 

1.84E+01 

-8.20E-01 

6.46E+01 

2.10E+01 

4. 60  E-01 

-2.76E-03 

2.36 

4.65 

Jun.  21 

9.9IE+0I 

3.63E+01 

-1  75E+00 

I.39E+02 

4.42E+0I 

9.02  E-01 

-5.44E-03 

2.17 

4.26 

Sep  21 

4.79E+01 

1.5  IE+01 

-8.86E-01 

6.03E+01 

2.41E+01 

1.1IE-01 

-6.65E-04 

1.65 

3.39 

Dec  21 

5.39E+00 

2.73E+00 

-1.22E-0I 

8.47E+00 

3.00E+00 

6.44E-02 

-3.85E-04 

2.48 

4.90 

pdf 

Mar  21 

3.46E  +  01 

1.06E+01 

-9.87E-01 

6.46E+01 

2.10E401 

1.49E-01 

-l. 00  E-03 

1.77 

3.45 

Jun  21 

6.96E+01 

1.85E+01 

-1.84E+00 

I.39E+02 

4.42E+01 

1.86  E-01 

-1.34  E-03 

1.61 

3.40 

Sep  21 

3.97E+01 

7.89E+00 

-8.65E-0I 

6.03E+01 

2.41E+01 

-1.50E-01 

8.26E-04 

1.24 

2.54 

Dec  21 

5.11E+00 

1.99E+00 

-1.70L-01 

8.47E+00 

3.00E+00 

3.85E-02 

-2.42E-04 

2.11 

4.12 

pJtr 

C<7 

Mar  21 

9.35E+00 

7.75  E+00 

1.67E-01 

-3.22E-13 

1.90E-12 

3.11  E-01 

-1.76E-03 

4.58 

9.79 

Jun  21 

2.94E+01 

1.78E+01 

8.64E-02 

-4.81E-13 

2.74E-12 

7.I6E-01 

-4.11  E-03 

4.11 

9.03 

Sep  21 

8.13E+00 

7.23E+00 

-2.08E-02 

-2.76E-13 

1.38E-12 

2.60E-01 

-1.49E-03 

4.84 

10.27 

Dec  21 

2.81  E-01 

7.38E-01 

4.75E-02 

-1.80E-14 

5.71  E- 1 4 

2.59E-02 

-1.43E-04 

5.98 

14.20 

Same  as  Table  3,  bul  for  mid-lalilude  atmosphere  wilh  urban  aerosols.  The  Table  shows  excellent  fil  of  Ihe 
regression  model  and  2.1%  mean  error  (all  36  entries).  The  worsl  accuracy  (~6%),  is  for  radiation  fields  E[ fr  and 


F£%/El  at  21  December.  For  all  other  radiation  fields,  the  mean  absolute  error  is  <5%  (mean  of  1.86%  for  all  34 
entries). 
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Tabic  5.  Mean  | error  for  Latitude/Datc  Scaling  of  a  Reference  Atmosphere 
to  Four  Dates  and  Four  Atmospheres _ _ 


Location  and  time  scaling  for  regression  model  xilantmk 
for  atmospheres  A,  B,  C,  and  D 

35°.Ator.2l'')- 

->  y(latituik\Ja\  ) 

Latitude  25 

Latitude  35 

Latitude  45 

Day  of  Year 

Mean  error\  (%) 

Mean  \error\  (%) 

Mean  error  (%) 

A 

B 

c 

D 

A 

B 

c 

D 

A 

B 

c 

D 

£r  /  Ft 

'  rm«x 

3.0 

3.0 

18.5 

18.6 

2.0 

2.2 

13.7 

13.6 

15.9 

16.4 

2.2 

2.4 

pdtf  ,  r-Jif 
*-•0  '  max 

2.9 

2.8 

18.2 

18.2 

1.8 

2.0 

13.5 

13.4 

14.4 

14.7 

l.l 

1.1 

Feb  1 

pdtf  ,  pT 

1  mix  1  max 

2.8 

3.0 

2.5 

2.8 

16 

1.2 

1.7 

1.3 

8.5 

6.4 

8.5 

6.6 

et 

1.5 

1.8 

23.4 

23.5 

12.2 

12.9 

31.7 

32.1 

24.0 

24.9 

40.3 

40.9 

Edif 

2.0 

1.9 

20.6 

20.5 

2.4 

2.3 

22.0 

21.6 

5.3 

4.6 

23.5 

22.8 

pdir 

2.5 

2.4 

21.8 

21.8 

18.8 

18.6 

36.9 

36.8 

46.5 

46.2 

58.3 

58.2 

ft  / ft 

c,a  1  'max 

0.9 

1.0 

0.9 

1.0 

3.7 

3.6 

3.7 

3.6 

10.7 

10.6 

10.7 

10.6 

pdif  j  pdtf 
*^a  1  max 

0.9 

1.0 

1.2 

1.3 

3.7 

3.7 

3.6 

3.5 

10.8 

10.7 

10.7 

10.7 

May  1 

pdtf  j  pT 

1  max  '  max 

1.0 

1.0 

1.0 

1.0 

5.2 

5.7 

5.3 

5.8 

12.5 

12.7 

12.6 

12.7 

E\ 

8.2 

9.0 

8.1 

8.8 

4.3 

4.8 

4.2 

4.7 

2.1 

2.0 

2.2 

2.0 

r-iiif 

2.8 

2.9 

3.2 

3.3 

3.5 

3.6 

3.9 

3.9 

5.0 

5.1 

5.2 

5.2 

r-dtr 

10.8 

10.6 

10.9 

10.7 

6.8 

6.7 

6.8 

6.7 

1.7 

1.7 

1.7 

1.7 

ft  /  ft 

'  rm»x 

0.3 

0.3 

0.3 

0.3 

5.1 

5.1 

5.1 

5.1 

12.3 

12.3 

12.3 

12.3 

pdtf  ,  pdtf 

C'a  'max 

0.3 

0.3 

0.3 

0.2 

5.2 

5.1 

5.0 

4.9 

12.5 

12.4 

12.3 

12.3 

Aug  1 

pdtf  /  pT 

'max  'max 

1.4 

2.0 

1.5 

2.0 

6.9 

7.5 

7.0 

7.5 

14.5 

14.8 

14.5 

14.8 

et 

7.0 

7.8 

6.8 

7.6 

3.5 

4.1 

3.4 

4.0 

2.4 

2.1 

2.4 

2.2 

pdtf 
c a 

4.3 

4.4 

4.8 

4.9 

5.3 

5.4 

5.7 

5.8 

7  0 

7.0 

7.2 

7.2 

pdtr 

Z'a 

9.9 

9.7 

9.9 

9.8 

6.6 

6.5 

6.6 

6.5 

2.3 

2.2 

2.2 

2.2 

Et  f  ft 

*^a  1  'max 

3.1 

3.1 

18.7 

18.7 

1.4 

1.5 

14.4 

14.4 

12.4 

12.7 

3.8 

3.6 

rdlf  j  pdtf 
*^a  rmax 

3.0 

2.9 

18.3 

18.4 

1.3 

1.5 

14  2 

14.2 

1 1.4 

11.7 

4.2 

4.1 

Nov  1 

pdtf  /  pT 
'  max  '  max 

2  9 

2.9 

2.6 

2.7 

1.3 

0.9 

1.4 

1.1 

7.7 

6.0 

7.8 

6.1 

El 

0.3 

0.4 

22.1 

22.1 

9.8 

10  4 

29.9 

30.2 

21.5 

22  4 

38.7 

392 

Ef 

1.0 

0.8 

21.4 

21.3 

1.9 

1.8 

22.5 

22.3 

4.9 

4.4 

24.3 

23.6 

r  dtr 

3.8 

3.8 

19.2 

19.2 

13.1 

12.9 

32.3 

32.2 

37.8 

37.5 

51.6 

51.5 

Each  of  Ihe  six  quanlilies  are  scaled  (eq  16)  from  a  reference  atmosphere  al  a  reference  dale  and  latitude  to  four 
atmospheres  (A,  B,  C,  and  D)  al  four  dales  (1  February,  1  May,  1  August,  and  1  November)  and  three  latitudes  (25, 
35,  and  45°).  The  error  is  computed  over  all  parameter  spaces (h, albedo, visibility) .  Atmosphere  A  is  U.S.  1976  standard 
with  rural  aerosols.  Atmosphere  B  is  U.S.  1976  standard  w  ith  urban  aerosols.  Atmosphere  C  is  mid-latilude-summer 
atmosphere  with  rural  aerosols.  Atmosphere  D  is  mid-latitude-summer  atmosphere  with  urban  aerosols.  1  February 
and  1  November  are  fall-w  inter  season.  1  May  and  1  August  are  spring-summer  season.  Reference  atmosphere  is 
same  as  in  A,  B,  C,  and  D  atmospheres  al  spring-summer.  The  reference  chile  and  location  are  21  March  al  latitude 
35°.  Average  of  all  288  entries  is  -10%,  where  the  error  is  <5%  for  1 30  entries. 
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